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A new LISA simulator (LISACode) is presented. Its ambition is to achieve a new degree 
of sophistication allowing to map, as closely as possible, the impact of the different sub- 
s ! ■ systems on the measurements. LISACode is not a detailed simulator at the engineering 



level but rather a tool whose purpose is to bridge the gap between the basic principles of 
LISA and a future, sophisticated end-to-end simulator. This is achieved by introducing, in 
a realistic manner, most of the ingredients that will influence LISA's sensitivity as well as 
the application of TDI combinations. Many user-defined parameters allow the code to study 
different configurations of LISA thus helping to finalize the definition of the detector. Another 



5^ ' important use of LISACode is in generating time series for data analysis developments. 
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The LISA mission aims at detecting gravitational waves (GW) from space, using 3 detectors 
located on the summits of a triangle whose center orbits around the Sun, following the Earth, 



at a distance of ~ 50 x 10 6 km. The distance between spacecraft is of the order of 5 x 10 6 km. 
Located inside each spacecraft, two (or one of two) proof masses are kept on a geodesic orbit by 
a drag-free system (DFS): a system that compensates, through the use of micro-thrusters, for the 
external forces acting on the spacecraft by keeping it at its proper position with respect to the 
proof masses. 

The triangular configuration, which depends on the orbits each spacecraft follows, is not rigid 
and will both rotate around its center and "flex" (i.e. the distance between spacecraft is time 
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dependent) with a period of one year. On a much smaller length scale, picometers, and at fre- 
quencies between 10 -4 and 10 -1 Hz (the current instrumental bandwidth), the distance between 
spacecraft will be modified by the passage of a gravitational wave (GW). The detection of such GW 
will be achieved by interferometric measurements giving the phase (frequency) difference between 
local and distant (i.e. coming from another spacecraft) laser beams on each spacecraft. There are 
therefore, 6 independent laser beams and 12 interferometric measurements: 6 between local lasers 
and 6 between local and distant lasers. A detailed description of LISA, of its principles and of the 
different noise sources, is given in the LISA Pre-Phase A Report [l]. 

The LISA detector system is complex and its sensitivity depends to a large extent on the 
different noise contributions coming from the DFS and the interferometric measurements. For 
these reasons and because a laboratory replica of the system is not totally achievable, the exact 
evaluation of the performances of LISA can only be studied with computer simulations of the 
different processes involved. Two such simulators have been elaborated in the US: Synthetic LISA 

and LISA Simulator [J. 

The present status of the LISA project warrants that precise evaluations of its performances 
be performed, compared and checked through different simulation softwares. The software imple- 
mentation of the LISA method and, for example, the introduction of the noise, the measurement 
of the phase, the proper time stamping of the measurements and the reduction of the laser phase 
noise by the TDI method are all non trivial and may depend on the manner in which different 
"technologies", not always perfectly defined at the present status of the project,are implemented 
at the software level. Of the two simulators mentioned above, Synthetic LISA 2( is the one that 
appears closest to mapping these technologies. However it remains necessary to develop similar 
software simulators in order to push this implementation further and to check that the quality and 
the sensitivity of the LISA detector is independent of the structure and numerics of these codes. 

This paper presents such a new LISA simulator, LISACode, whose ambition is twofold: i) a 
modular, semi-technical, implementation of the different ingredients, and particularly the different 
noise sources and the phasemeter, which will affect the quality of the measurements and ii) the 
application of the TDI (Time Delay Interferometry) , on the phasemeter outputs, which allows to 
extract the GW signal from the overwhelming laser noise. The main inputs of the code are user- 
defined or built-in GW sources and the various parameters which define the noise levels and the 
orbits of the spacecraft. The main outputs of LISACode are time sequences, with a user-defined 
sampling step, of the basic (phasemeter) measurements and/or time sequences of user-defined TDI 
generators. 



After this introduction, the second section of this paper presents the structure and different ele- 
ments of LISACode. The third section gives some examples of studies that LISACode can be used 
for and a last section presents a number of GW pre-defined in LISACode. Two appendices describe 
the reference frames used in LISACode and the CPU performance of the code in various configura- 
tions. The LISACode software package and a more detailed description, can be downloaded from 
the LISA-France web site at http://www.apc.univ-paris7.fr/LISA-France/analyse.phtml. 

II. LISACODE 

LISACode is not a detailed simulator at the engineering level but rather a tool whose purpose 
is to bridge the gap between the basic principles of LISA and a, future, sophisticated end-to-end 
simulator at mission level. This is achieved by introducing, in a realistic manner, most of the 
ingredients that will influence its sensitivity. The end product of LISACode are the time series of 
the phasemeter outputs and their fluctuations arising both from the instrumental noises and from 
the statistical fluctuations inherent to any measurement process. 

A. The Structure of LISACode 

The schematic structure of LISACode is presented in figure HJ The core of the program which 
simulates LISA is shown inside the box. The GW calculations and the application of TDI, although 
they are part of the code, are not part of LISA and can be considered as external modules. The 
basic output of the code are time series of the phasemeters. These time series can be calculated 
with a variety of different noise inputs with or without the presence of a GW. These time series 
can be post-treated by a number of pre-defined or user-defined TDI combinations. At the end of 
this section, various sensitivity curves for different TDI combinations will be presented showing 
the numerical accuracy that can be obtained by LISACode. 

As stated above, the aim of LISACode is to follow, as closely as possible, the technological 
design of the future detector. There are, at least, two ingredients that are fundamental in this 
respect : the different noises that intervene and limit the sensitivity of LISA and the precision of 
the measurement of the phase (or frequency) by the phasemeter module. Associated to this are 
the USO (Ultra Stable Oscillator, the reference clock on each satellite) which will impact both on 
the time-stamping of the time series and on the precision of the phasemeter. We discuss below 
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FIG. 1: Schematic diagram of the LISACode software simulator. 



some of the more important features of LISACode. 



B. The Orbits and photon flight paths 

The orbits of the three spacecraft are calculated using the analytical formulation of 4|]. These 
orbits include both the rotation of the LISA constellation and the flexing of the distance between 
the spacecraft. At the origin of time (t = 0), the center of gravity of the LISA constellation points 
towards the vernal point. The photon flightpath between two satellites can be computed at several 
orders in an expansion in powers of the Sun's gravitational potential, giving increasing degrees 
of accuracy Qj. At the lowest order, the light propagation time is computed using the distance 
between the spacecraft. Further precision can be obtained by taking into account the Sagnac effect, 
the aberration and other General Relativity effects Q]. 



The Noises 



In LISA, the noises can be divided in two groups depending on wether they are transported 
from one satellite to another (allowing their suppression, to a certain level, by the application of 
TDI) or wether they impact locally (and cannot be suppressed). 



• External noises: laser noises that are transported from one satellite to another and which 
will be suppressed by TDI. At this moment in time, the laser noise is considered to be a 
bandwidth limited white noise whose power spectral density is 30Hz/v / Hz 6j. It is almost 
certain that the effective laser noise will present a more complicated frequency dependence, 



notably at low frequencies. This will be implemented as soon as more information on this 
subject is available. 

• Internal noises: inertial mass, Ultra Stable Oscillators (USO), shot and optical path noise. 
These will be transmitted through the phasemeter transfer function and will impact on the 
precision of the measurements. 

Table U gives the different inputs to the internal noise for a phasemeter. The values of the second 
column are those given in reference The third column indicates the noise used, by default, in 
LISACode. These values can be modified, for each individual phasemeter, by the user. Note that 
the noise in lines 3 to 8 are, in LISACode, summed up quadratically, for each phasemeter, into a 
unique term. The last column of table U gives the assumed dependency of the different noise with 
the laser power and the distance between satellites. The residual laser noise does not appear in 
LISACode as this was assumed to represent the residual noise after suppression of the laser noise 
by TDI Q] which is treated explicitly by LISACode (see section ITTF|) . 

When translated to fractional frequency fluctuation (Ais/v unit) all but the last noise given in 
this table have an / dependence at the phasemeter level. In the same units, the inertial mass noise 
has an 1/f dependence. 

Q 

A number of remarks have to be made on this table. Since it was issued, ESA [7( has made a 
more refined study of the allowed noises and has subdivided them into different components. At 
the phasemeter level, all the measurement noises (using the terminology of TJ) have a sj f 2 + a/ f 2 
frequency dependence, whereas the acceleration noise has a y/b/ f 2 + f 2 /c + a// 6 frequency de- 
pendence. This subdivision of different noises and their frequency dependence are being imple- 
mented in LISACode but the results presented in this paper maintain the definition and frequency 
dependence given in table E 



D. The Phasemeter 



In LISA, measurements are performed at the phasemeter level (see fig. [T]). A phasemeter sums 
up the effects of the gravitational strain and of the various noises. There are two phasemeters for 
each optical bench (see fig. [2]): one (sj) measures the phase between the external laser and the 
local one, the other (n) measures the phase between the local laser and the local laser of the other 
optical bench in the same spacecraft. For example, Equations [U El E] and H] give the structure of 
the signals for the four phasemeters of spacecraft 1: 



TABLE I: Noise allocation budget based on Table 4.1 of The second column gives errors listed in the 
PrePhaseAReportjl]]. The third column gives values used in LISACode in the same unit as the second 
column. The fourth column gives the default errors included in LISACode and their frequency dependence. 
These values can be modified by the user. 
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where the unprimed (resp. primed) quantities refer to the link i— ► i+1 (resp. i— > i-l)(see figureE]). 
For any function of time /, Dif = f(t — ti), D'J = f(t — t^), are the time delay operators, and the 
delays tj,t£, the photon flight times in the two possible directions. Di takes into consideration the 
travel time in the direction 3 — > 2 — > 1 and D\ in the direction 1 — > 2 — > 3. sf w and s'f W are the 
gravitational signals on the laser link respectively from the spacecraft i + 1 to the spacecraft i and 
from the spacecraft i — 1 to the spacecraft i. s Meas - nmse groups the shot noise and other optical 
path noises (see table H]) received on phasemeter Sj. p l . aser - nmse are the laser noise and §A cc - nmse j s 
the inertial mass noise (see table [I]). All these signals are measured in relative frequency fluctuation 
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unit (Au/u unit). 




The impact of the phasemeter measurements and errors on the final quality of the time series is 
a sensitive issue and it is a specific feature of LISACode that it attempts to tackle this technological 
problem as precisely as possible in view of our present knowledge the phasemeter. Traditionally, 
the phasemeter outputs will be transmitted to earth at a rate which may vary from 1Hz to 3Hz. 
In order to eliminate the high frequency component (i.e. frequencies above 0.5 the transmitted 
data rate) the phasemeter must be associated with a filter. In LISACode this filter is implemented 
in the following way: once the sampling rate of the phasemeter output is chosen by the user, the 
time series are calculated at a multiple (user-defined) of this frequency and a filter is applied to 
suppress the frequencies above a tenth of the sampling rate. The outgoing signal will therefore 
still have a (small) residual high frequency component (i.e. above the cut-off frequency) which will 
impact on, for example, the application of TDI. Although the final filtering properties associated 
to the phasemeter may depart from this implementation, it is considered that such features should 
be present in this kind of simulation in order to evaluate the impact it may have on the final 
performance of the detector and of the analysis. Figure [3] of section III Fl shows (red curve above) 
an example of the phasemeter +filter output for a sampling rate of 1Hz and illustrates the cut-off 
obtained above 0.1Hz . 



E. The Reference Clocks 



The USO (or reference clock) noise impacts the measurements at two levels. As the USO (one 
in each satellite) will be used for the phasemeter module by giving it a reference signal, its noise 



will impact on the phasemeter output. This is accounted for by the corresponding noise given on 
the fourth line of Table [H The USO is also used to time — stamp the time series and will thus 
introduce an error in the application of TDI. In LISACode, the time stamping is affected both by 
this noise as well as by a, user-defined, offset and a drift with respect to barycentric time. At the 
phasemeter level, only the USO noise is considered. The more complicated effect of a drift is under 
study. 



F. TDI and the LISA Sensitivity Curve 

Due to the large distance between the spacecraft, the interferometric measurements are per- 
formed via independent (local and distant) lasers. The laser noise level is typically of the order of 
10~ 13 (in Av/v unit), whereas, the gravitational wave signal is closer to 10~ 21 . This implies that 
the laser noise must be suppressed by more than 8 orders of magnitude. Today, this is achieved 
by a numerical method called Time Delay Interferometry (TDI). An extensive description of 
the TDI method, and of the different combinations, can be found for example in [§[[1], 
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TDI is fully implemented in LISACode. 

The general philosophy of TDI is to linearly combine the phasemeter measurements obtained at 
different times and spacecraft in such a way as to eliminate numerically the laser noise. There are 
various linear combinations that allow for this. As an example the Michelson X (1 st generation) 
generator corresponds to the combination: 



X = -81- D 3 s' 2 - D 3 Dy s'i - D 3 DyD 2 , s 3 

(5) 

+s'i + D 2 < s 3 + D 2 >D 2 si + D 2 ,D 2 D Z s 2 > 

where Sj is a phasemeter output (see eq. [I]) and D{ is a delay operator based on the time of flight 
of the the laser over the link Lj: 

D iX (t) = x(t - Life) (6) 

Such combinations will tend to zero with a residual value which will depend on the level of the 
experimental noises (see section III CP and on the accuracy with which the phasemeter time series 
can be interpolated at the corresponding delay times. 

LISACode includes a number of predefined TDI combinations, for a fixed (1 st generation) and 
a moving (2 nd generation) LISA: X, Y ,Z , a, (3, 7 and £ as well as Relay, Monitor and Beacon. 



A numerical application of TDI by LISACode is illustrated in figure [3] which shows the Power 
Spectral Density (PSD) of the phasemeter time-series signal (red curve above) before TDI post- 
processing, and the PSD of the Michelson TDI X 2 nd generation (black curve below). The GW 
peak is clearly visible at of 10 -3 Hz. The cut off observed at 0.1 Hz is due to the phasemeter's 
filter as discussed in Section III PL The inset in fig. [3] shows the structure of this peak that reflects 
the doppler modulation of the initial frequency due to the movements of LISA over a year. 
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FIG. 3: Power spectral density of the phasemeter si (red curve above) and of the Michelson TDI X 2 nd 
generation (black curve below). The GW is at a frequency of 10~ 3 Hz and has an amplitude of 10~ 19 . The 
inset shows an expansion around 1 mHz (Color on line). 

A way to test the numerical accuracy of LISACode is to reproduce the sensitivity curve of LISA 
as obtained for a given TDI combination. The sensitivity curve represents LISA's sensitivity, as 
a function of frequency, for an isotropic distribution of sources with random polarization. The 
sensitivity s| corresponds to the square root of the ratio of the noise spectral density PSD Noise to 
the duration of the measurement (1 Year) divided by the gravitational wave TDI response Resp^fiy, 
and for a signal to noise ratio of 5, i.e.: 



PSD Noise L 

V lYear X Resp r ^ [ ' 

This function is traditionally calculated using X TDI 1 st generation, for a rigid and non-rotating 



LISA, for which a semi-analytical calculation is possible [8l. 1 121.1 131]. This semi-analytical calculation 
is performed without laser noise, assuming that it is exactly cancelled by the TDI preprocessing. 

As calculating with LISACode the sensitivity curve for a precise isotropic distribution of ran- 
domly polarized sources would imply a prohibitive CPU time, we have first determined a limited 
set of sources which yield a sensitivity curve very close to the one given by the exact calculation: 
an agreement better than 10% for frequencies below 0.1Hz. The limited set is defined by 44 
(incoherent) sources whose parameters can be obtained fromfl^ 

Figure H] shows the exact semi-analytical (solid red line) calculation and the numerical calcu- 
lations performed with LISACode, including laser noise, as (black) dots. As can be seen, the 
sensitivity obtained from the limited set reproduces very accurately the exact semi-analytical re- 
sults. These calculations can be compared to results presented, for example, in and [§| (a Web 
121 ] site also allows to calculate similar sensitivity curves). These results allow us to conclude that 
this subset of sources can be used to test sensitivity calculations on a wide frequency range and 
that the implementation of TDI in LISACode properly eliminates the laser noise (see also figure [6] 
for an illustration of this). In all the following calculations, this subset of sources will be used to 
generate the different sensitivity curves. 
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FIG. 4: The comparison between the analytic sensitivity curve (solid red line) for an isotropic distribution 
and the sensitivity curve obtained with LISACode (black dots) for the same set, see section Hi Fl and ref [3] 
(Color on line). 
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FIG. 5: The comparison between the sensitivity curves obtained for a moving LISA. The (red) crosses show 
the sensitivity curve using X TDI I s * generation and the (black) dots the sensitivity curve using X TDI 2 nd 
generation. Curves calculated with a set of 44 sources, see section Hi Fl and ref (Color on line). 

Figure compares the sensitivity curves obtained for a moving LISA. The (red) crosses show 
the sensitivity curve using X TDI I s * generation and the (black) dots the sensitivity curve using 
X TDI 2 nd generation . The impact of the motion of LISA, and the imperfect annulation of the 
laser noise for the I st generation TDI is clearly observed. The application of the 2 nd generation 
combination corrects this and recovers a sensitivity curve very close to the one shown in figure HI 

Figure [6] shows the noise budget of the X TDI 2 nd generation sensitivity curve. At low fre- 
quencies, the accelerometer noise (i.e. inertial mass noise, c.f. Table [J) dominates whereas at high 
frequencies the measurement noise, which includes the shot noise, is the most important factor. As 
can be seen, the application of TDI, (assuming a laser with a power spectral density of 30Hz/VHz), 
puts the residual laser noise below the other contributions (see Section III Cj) . Note that this result 
depends on the precision with which the interpolation of the phasemeter time series is performed 
when using the delay operator (c.f. eq. [6]). In order to cover the whole frequency range we 
have concluded, in agreement with M.Vallisneri , that a 20* /l order Lagrange interpolation is 

necessary. 
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FIG. 6: The budget of each noise contribution to the sensitivity curve of TDI X 2 nd generation. The optical 
path noise including the shot noise are shown by (red) squares, the inertial mass noise by (green) crosses, 
the residual laser noise by (blue) plusses and the sum of all by (black) dots. Curves calculated with a set of 
44 sources, see section Til F I and ref Q (Color on line). 

III. APPLICATIONS OF LISACODE AND PRE-DEFINED GW 

The LISA community is in the process of finalizing the design of the detector and there are 
numerous areas in which a simulation software, like LISACode, can help to define the optimal 
configuration and decide between different options. In the following, we apply LISACode to some of 
these questions: the effect of changing the nominal LISA armlength, the accuracy in the knowledge 
of the flightpath time on the sensitivity of LISA and the TDI combinations that can be used in 
the case of a missing link. 

A. Modifying the Armlength: from 5 x 10 9 m to 2 x 10 9 m 

The scale of LISA (nominal armlength) has been a subject of debate within the LISA commu- 
nity as it has sometimes been advocated that a smaller LISA could result in a reduction of its 
financial budget. LISACode can address the scientific aspect of this question which is relative to the 
dependence of the sensitivity curve on the nominal armlength. Figure [7] compares the sensitivity 
curves for armlengths of 2 x 10 9 m and 5 x 10 9 m. One observes that the sensitivity for the shorter 



armlength is shifted towards higher frequencies and that the "oscillation" peaks, conditioned by 
the armlength value, are positioned at the corresponding frequencies. 
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FIG. 7: Sensitivity curves for two armlengths: (black) dots for a 5 x 10 9 m armlength and (red) plusses 
for a value of 2 x 10 9 m. Calculation are performed in the framework of TDI X 2 nd generation. Curves 
calculated with a set of 44 sources, see section IPTF! and ref Q (Color on line). 



B. Inaccuracy on the delays 

Application of TDI requires a very good precision on the travel time of the lasers beams between 
spacecraft (see Equation \5\ and the discussion in This travel time will be known with a given 

precision and the sensitivity of LISA will depend on this precision. In order to test the sensitivity 
as a function of the precision, an offset has been added to the exact delay used in a given TDI 
combination: 

Dtdi = D exact + AD Q ff set (8) 

This offset, AD Q ff se t, is assumed to have a magnitude of about 1 //sec, i.e. a photon flight path of 
about 300 m. Figure [8] shows the evolution of the sensitivity for TDI X 2 nd generation for an offset 
of 0, 0.5 and 1 fisec respectively. As the level of the laser noise is of the same order of magnitude 
as the level of the other noises, the sensitivity curve goes up as a function of the offset due to the 
increasing error in the interpolation scheme produced by the inaccuracy of the delay. 




FIG. 8: The sensitivity curves for different offsets on the laser travel time used in TDI: (black dots), 
0.5 (red plusses) and I fisec (blue crosses). Calculations are performed for TDI X 2 nd generation. Curves 
calculated with a set of 44 sources, see section Hi Fl and ref [14[ (Color on line). 



C. The effect of a missing link 

A question that the LISA mission may face is that of the loss of a link (from spacecraft i to 
j, but still keeping the link of j to £) between spacecraft since the question of the reliability of 
the links has also an important impact on its financial budget. Several TDI combinations have 
the property of allowing for the loss of such a link. Figure [9] shows the Beacon, Monitor and 
Relay combinations 



111 ] and the Michelson 1 st generation which requires all the links between one 
spacecraft and the two other. If all the links are available, all three Michelson (A, Y and Z) would 
be conserved whereas, in the case of the loss of a link , one would conserve only one Michelson and 
the Beacon, Monitor and Relay combinations (Note that the Beacon and Monitor are numerically 
almost indistinguishable). The loss of precision in determining the physical parameters of a source 
when a link is loss has still to be determined. 



IV. PRE-DEFINED GW IN LISACODE 



The LISACode software package includes a number of predefined types of gravitational waves 
which are rapidly described below. 
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FIG. 9: The comparison of sensitivity curves for TDI combinations Beacon (red squares), Monitor (blue 
crosses) and Relay (green plusses) with X 1 st (black dots). Curves calculated with a set of 44 sources, see 
section III Fl and ref 14[ (color on line) 



A monochromatic wave: 



This type of wave is essentially a "mathematical" object which is directly defined in the Canon- 
ical Reference Frame (see fig. [T2l) . Therefore the inputs to LISACode will be the amplitudes /io+ 
and ho* of the two polarization components, the frequency / and 0o+ and 4>q x the initial phases 
of the two polarization components. The conventions for the strain calculation are: 

h C RF+{t) = /io+sm(2vr/t + ^ +) 

(9) 

hcRFx{t) = h 0x sin (2tt ft + (p 0x ) 



B. GW signal from a binary with fixed frequency: 

A binary of fixed frequency goes one step further. LISACode defines the amplitude by the 
physical parameters of the system in the canonical reference frame (see Appendix [X] and fig. fl~2j) : 

= A (l + cos 2 i) cos (Ait f orb t + 4> ) 
= 2 A cos i sin (47r/ or fet + <po) 




with 



mtot = m 1 +m 2 (11) 

r= &?T (i2 » 



where ra\ and mi are the masses of the two objects, f or b is the orbital frequency, i the inclination, 
00 the initial phase and r the distance from the source to LISA (c is the light velocity). 



C. GW signal from a coalescent binary, computed in the Post-Newtonian approximation: 

The input parameters are the same as above except for the additional parameter T coa i escence 
which defines the coalescence time. Figure [H)] shows an example (time-frequency matrix for a 5 
month period) of such a calculation superposed on the LISA instrumental noise. Reference [3] 
gives an explanation of the model which is used. 




Time(s) jcq" 

FIG. 10: Time frequency figure of the Michelson TDI (X 2 nd generation) for a binary system computed 
in Post-Newtonian approximation at 2.5 PN during 5 months (10 6 - 10 6 solar masses, coalescence time 4 
months, inclination 90 degrees and distance 100 Mpc). The color scale represents the log of the amplitude 
(Color on line). 



D. Input of a GW from a Text file. 



In the pre-defined GWs described above, the direction of the angular momentum of the emitting 
system does not depend on time and, therefore, neither do the unit vectors p and q and, hence, 
the inclination and polarization angles. 

In more complicated cases, such as for EMRIs (Extreme Mass Ratio Inspiral) 1?J , this is not 
the case and therefore besides hcRF,+ and hcRF,x, the inclination and polarization angles would 
have to be given as a function of time. To simplify the inputs, LISACode requests that the strain 
be expressed directly in the "Observational Reference Frame" (ORF, see fig. [T2|) . The requested 
(ASCII) input file should therefore contain, for each time step, a line giving the time, hoRF+{t) 
and hoRFx(t). The data can be given with any constant time step as LISACode will linearly 
interpolate this data to obtain data sampled at the correct, user-defined, phasemeter time step. 
More details on the input file can be found in Figure [TT] shows an example (Time- Frequency 
matrix for a one year period, including the detector noise) of harmonics produced by an EMRI 
with a mass ratio of 10 5 , a spin ratio of S/M = 0.8, at a distance of 100 Mpc. 



,x 10 




Time (s) 



FIG. 11: Time frequency figure of the Michelson TDI {X 2 nd generation) for EMRI during 1 year (10 - 10 6 
solar masses, black hole spin ratio 0.8 and distance 100 Mpc ). The color scale represents the log of the 
amplitude (Color on line). 



V. CONCLUSION AND PERSPECTIVES 



We have presented a new simulator (LISACode) whose ambition is to achieve, through a very 
modular structure, a new degree of sophistication allowing to map, as closely as possible, the 
precise impact of the different sub-systems on the precision of LISA measurements. LISACode 
is considered as a tool whose purpose is to bridge the gap between the basic principles of LISA 
and a, future, sophisticated end-to-end simulator. This is achieved by introducing, in a realistic 
manner, most of the ingredients that will influence its sensitivity. Most of the parameters of a 
LISA configuration (orbits, noise functions,...) can be user-defined. LISACode can use number of 
pre-defined gravitational waveforms but more complicated waveforms can also be introduced via 
text files. As the suppression of the laser noise is a very important aspect of LISA, LISACode allows 
to implement a number of pre-defined TDI combinations and allows the user to construct more 
specific combinations if needed. In its most complete form, LISACode is relatively time consuming 
but different and faster versions of LISACode can be constructed and used, for example, for data 
analysis (see Appendix [B]) . 



APPENDIX A: REFERENCE FRAMES AND THE RESPONSE OF LISA TO A GW 



There are a number of reference frames which have to be defined in order to properly understand 
and generate the response of LISA (see for instance : 

m 

The first one is the Barycentric Ecliptic Frame (BEF), centered on the Sun. Its X-axis is 
defined by the direction pointing towards the vernal point. The Z-axis is defined by the direction 
of the ecliptic pole and the Y-axis completes the reference frame. 

In this frame, the source direction n is defined by the declination (3 and the longitude A, as 
shown in fig. fl2l This is the reference frame that is used by LISACode for computation. 

Using this direction, an orthonormal frame (6, <j), n), called the Observational Reference 
Frame (ORF) is defined by: 

^cos (3 cos a\ / sin (3 cos a\ / sin A ^ 

c- dn 

sin (3 sin A , 1 
i — cos (3 J \ u j 
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cos (3 sin A 
i sin (3 J 
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— 1 dn 

cos (3 d\ 



cos A 
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The 6, <p plane is called the observation plane. Also defined by the source direction, the Canonical 
Reference Frame (CRF: p, q, n) is defined. The unit vector p is defined as the direction of the 




y : verrmi point 

FIG. 12: Schematic description of GW parameters and reference frames. The direction of the source is 
located by the declination (3 and the ecliptic longitude A. (6, <f>, n) define the Observational Reference 
Frame , i.e. the frame constructed from the direction of the source. In this frame 6 is on the meridian, (p, q, 
n) define the Canonical Reference Frame, ip is t ne polarization angle, i. e. the angle between the Canonical 
Reference Frame and the Observational Reference Frame (Color on line). 

major axis of the ellipse obtained by the projection of the orbit on the observation plane. The 
definition of q follows from the product p x n. 

Geometrically, and in the case of a circular orbit, the unit vectors p and q (labelled a and b in 
[{]]) are aligned with the major and minor axis of the projection of the orbit on the 9, 4> plan. In 
terms of polarization, p and q define the reference frame within which the gravitational strain is 
given by: 



(h 



hij(t) 



CRF+(t) hcRFx(t) 
hcRFx(t) -h C RF+(t) 

0/ 



(A2) 



Where hcRF+(t) arid hcRFx (t) are the time evolution of the plus and cross polarization components 
in the canonical frame (CRF). 

The inclination angle i is defined as the angle between L (orbital angular momentum of the 
emitting system) and the source direction n. 



The polarization angle, ip, is defined as the angle between 9 and p, i.e.: 

p = cos ib — sin ib 6 

^ ^ 1 (A3) 

q = sin ib + cos ib (ft 

It is important to notice that in many cases, p, q, the polarization angle and the inclination 
angle will be time dependent and will therefore imply special care in defining the gravitational 
wave strain for LISACode. An example of this is given in section IIVDI for an EMRI where a GW 
is loaded to LISACode via a text file. 

The strain of a GW produces a phase fluctuation in the laser link between two spacecraft. Since 
this phase fluctuation is related to the frequency fluctuation, the LISA response involves the six 
frequency fluctuations of the laser links. The relative GW frequency fluctuation for a laser received 



by spacecraft B from spacecraft A, is directly connected to the gravitational strain 

Vopt 



0] 



by: 



/.x _ H AB (t + n-rt) - H AB (t + n-rX- L A b/c) . . 

AB 2(1 + n-nA~B) [ ' 



with 

H A B(t) = h RF+{t) £,+,AB + h RFx (t) £,x,AB (A5) 

where n is the unit vector of the source's direction, r~X and Tb are the position of the spacecraft 
in the barycentric reference frame, L AB /c is the time of propagation of a photon between the two 
spacecraft, raJe is the unit vector along arm AB, hoRF+(t) and hoRFx(t) the time evolutions of 
the GWs two components in the observational frame (6, (ft, n). £,+ iAB and (, X ,AB are the two 
antenna pattern functions: 



£ xi = 2 \J3 ■ Hi ) [X-n 

The relation between the GWs two components in the observational frame (ORF: 8, (ft, n) and 
the components in the canonical frame (CRF: p, q, n) is: 

hoRF+ = cos(2V0 h C RF+ + sin(2V>) h C RFx 

(A7) 

hoRFx = -sin(2V>) h C RF+ + cos(2ift) Ii C rfx 
APPENDIX B: CPU CONSIDERATIONS 



In its complete configuration, LISACode aims at generating time sequences of data with realistic 
noise functions and testing the impact of different LISA configurations: orbits, TDI combinations, 
etc... 



It soon appeared necessary to use LISACode within a data analysis strategy. For this, the 
CPU execution time is a critical issue sometimes overriding the priority of numerical precision. For 
instance, assuming exact cancellation of laser noises saves a significant amount of CPU time. Other 
options as interpolation of the orbits and using less efficient but simpler TDI combinations can also 
speed up the execution. Table HT1 gives some execution times, measured on a G5 biprocessor Mac 
(2.7 GHz) running Unix. The last line corresponds to using a simplified version, LISALight, which 
is compiled with modules extracted from the LISACode libraries in order to directly calculate the 
TDI combinations. Calculations are performed for a monochromatic wave of frequency lmHz. 

TABLE II: Execution time for LISACode simulations for one year with a time sampling of At sampling- 
Order 2 orbits take into account flexing, Sagnac, aberration and relativistic gravitational effects. Order 
orbits include flexing only. LISALight corresponds to a program compiled from the LISACode libraries but 
simplified as much as possible. Calculation were performed on a G5 biprocessor Mac (2.7 GHz, 5 Go RAM) 
running Unix. 
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